Condensed Fraction of an Atomic Bose Gas Induced by Critical Correlations 
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We study the condensed fraction of a harmonically-trapped atomic Bose gas at the critical point predicted by 
mean-field (MF) theory. The non-zero condensed fraction /o is induced by critical correlations which increase 
the transition temperature Tc above . Unlike the Tc shift in a trapped gas, /o is sensitive only to the 
critical behaviour in the quasi-uniform part of the cloud near the trap centre. To leading order in the interaction 
parameter a/Xo, where a is the s-wave scattering length and Aq the thermal wavelength, we expect a universal 
scaling /o oc (a/Ao)^. We experimentally verify this scaling using a Feshbach resonance to tune a/Ao. Further, 
using the local density approximation, we compare our measurements with the universal result obtained from 
Monte-Carlo simulations for a uniform system, and find excellent quantitative agreement. 

PACS numbers: 03.75.Hh, 67.85.-d 



Some of the most interesting fundamental problems of 
many-body physics involve strong inter-particle correlations, 
and cannot be addressed by mean-field (MF) theories. Har- 
monically trapped ultracold atomic gases are promising can- 
didates for highly controllable "quantum simulation" of such 
intricate many -body scenarios However, for testing the 
existing theories of spatially uniform systems, it is often im- 
portant to experimentally extract information on local proper- 
ties of a non-uniform trapped gas (see e.g. 

The effect of interactions on Bose-Einstein condensation 
of a dilute gas is a classic example of a difficult beyond-MF 
problem, which has challenged theorists for decades lIS HTTl . 
It is also an example of a situation where harmonic confine- 
ment both quantitatively and qualitatively alters the physics 
lITStiZTl . For a uniform gas the interaction shift of the crit- 
ical temperature Tc cannot be calculated to any order in the 
interaction strength using perturbation theory, owing to strong 
correlations that develop near the critical point. On the other 
hand, non-uniformity of a trapped atomic gas results in a sig- 
nificant MF shift of Tc ifTSl . More importantly, it diminishes 
the more interesting beyond-MF effects, in essence because 
near Tc only a small fraction of the cloud is actually in the 
critical regime (see Fig. [T}. Only recently have the beyond- 
MF effects on condensation of an atomic gas become experi- 
mentally accessible |27|. Many questions remain open since 
beyond-MF effects in a uniform and a trapped gas have differ- 
ent dependence on the strength of interactions, and quantita- 
tive connections between the two are highly non-trivial. 

In this Letter, we study the condensed fraction (/o) of an 
atomic Bose gas at the critical point predicted by MF theory. 
By definition /o vanishes within MF theory, and directly mea- 
sures the effect of critical correlations which shift Tc above 
T^^. Moreover, while the Tc shift itself strongly depends on 
the global properties of a non-uniform gas, fo measurements 
directly probe the quasi-uniform critical region near the cen- 
tre of the trap. To leading order in the strength of interactions 
we predict a universal scaling fo cx (a/Ao)'', where a > is 
the s-wave scattering length and Xq the thermal wavelength at 



the ideal gas critical temperature T^. Using a Feshbach res- 
onance in a ■^^K gas to tune q/Xq, and accurately measuring 
condensed fractions in the range 0.1 — 1%, we experimentally 
verify this prediction. Further, we directly relate our measure- 
ments to the universal critical behaviour seen in the classical- 
field Monte-Carlo simulations of a uniform system ||28l , and 
find excellent quantitative agreement. 

In Fig.[TJa) we illustrate the difference between the beyond- 
MF shifts of the critical point in a uniform and a trapped sys- 
tem, and in Fig. [TJb) the expected scaling of the condensed 
fraction at the MF critical point. For visual clarity, here we fix 
the temperature of the gas and consider the interaction shift of 
the critical density Uc (in the centre of the trap) and the crit- 
ical atom number Nc- Surprisingly, the beyond-MF shifts of 
Uc and Nc are not directly related to each other The quadratic 
beyond-MF Nc shift is of direct relevance to the experimen- 
tally pertinent case of a trapped gas, but from the point of view 




FIG. 1. (color online) Beyond-mean-field effects near the critical 
point in a harmonically trapped Bose gas. (a) For a fixed tempera- 
ture, the density distribution at the critical point N = Nc < N^^ 
(solid blue line) is compared with the mean-field prediction (dashed 
red line). In the trap centre we expect n^^ — ric oc a/ Ao, characteris- 
tic of the critical behaviour in a uniform system. However the experi- 
mentally measured A'^c shift, A'^^^ — A^'c oc ( a/ Ao)^, is dominated by 
the density shift outside the central critical region, and is not directly 
related to the Uc shift, (b) If is increased to A''^^ > A^c, a small 
condensate induced by critical correlations forms within the critical 
region of size oc a/Xg. The condensed atom number A^o oc (a/Ao)* 
directly relates to the critical density shift Arte oc a/Xo. 
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of the theory of critical behaviour the hnear ric shift is actu- 
ally more interesting. Here we will show how to use a trapped 
atomic cloud to experimentally obtain information about the 
critical behaviour in a uniform system. 

We first outline some general scaling arguments, then 
present our experimental results, and finally return to a quanti- 
tative comparison of our measurements with the theory based 
on the classical-field Monte-Carlo simulations of Ref. |28 1 for 
a uniform system. 

In a uniform system, ideal-gas condensation occurs at a 
chemical potential /i^ — 0, and a critical phase space den- 
sity = C(3/2) ~ 2.612, where ( is the Riemann func- 
tion. In an interacting gas there is no Tc shift at MF level, 
i.e. T^^ = T^. To leading order in a/Xo < 1 the expected 
beyond-MF shift is given by HTHTtI 



trapped gas fZTl qualitatively mirrors Eq. ([2]): 



at; 
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a 
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where AT^ = - T°, and c w 1.8 (TTirni. Equivalently, 
the shift at constant T is Anjn° « -(3/2)Arc/Tc^. 

An important point is that, at both MF and beyond-MF 
level, the interactions differently affect (or equivalently n.c) 
and the critical chemical potential jic- The simple MF shift 
/3^MF ^ 4(^(3/2)a/Ao, where /3 = l/fceT, has no effect on 
condensation, and to lowest beyond-MF order |29|: 



Ao 



(2) 



The qualitative difference between Eqs. ([T]) and (|2]i highlights 
the fact that the problem of the shift is non-perturbative and 
near criticality the equation of state does not have a regular 
expansion in fi (otherwise one would get Aric oc /i*^^ — Hc)- 

In a harmonically trapped gas is defined for a given atom 
number TV, rather than for a given density n. For an ideal gas 
fcBTP = hoj [iV/C(3)]^/^ where C(3) « 1.202. Within the 
local density approximation (LDA) one expects the uniform- 
system results for jic and jic to apply in the centre of the trap, 
r = 0. Elsewhere in the trap the local chemical potential is 
/i(r) = /x(0) — mcLi^r^/2, where m is the atom mass and w the 
trapping frequency. The result for the Tc shift however does 
not carry over so easily to the trapped case; the basic reason 
for this is that at Tc only a small fraction of the non-uniform 
cloud is actually in the critical regime. The size of the central 
critical region is Tc ^ {a/Xa)RT, where Rt = ^/ksT/ 
is the thermal radius |22|. Combining this with Auc ^ a/Ao 
implied by Eq. ([TJ, we obtain a very small beyond-MF shift of 
the critical number of atoms within the critical region, of the 
order (a/Ao)^. However the interaction shift of Hc affects the 
density everywhere in the trap. One famous consequence of 
this is the negative MF shift of Tc in a harmonically trapped 
gas ifTSi : while = n", repulsive interactions broaden the 
density distribution so that N^'^ > N°. 

More generally, the experimentally observed Tc shift in a 
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Here bi « —3.426 is an analytical, strictly MF result ifTSl . 
and 62 = 46 ± 5 was measured in ll27]l . The quadratic Tc 
shift depends on the beyond-MF correlations, but does not di- 
rectly correspond to the lowest-order (beyond-MF) Tc shift in 
a uniform system [Eq. ([T]i], which should enter the in-trap Tc 
only at the (a/Ao)* level. We can qualitatively understand the 
similarity of Eqs. (|2]i and ([3]l by noting that; (i) Away from 
the critical point the equation of state is regular in /i and the 
local density shift is simply proportional to ^c, at both MF 
and beyond-MF level ||30]| . and (ii) The contribution to Nc 
from the non-critical region outweighs the contribution from 
within the critical region by a large factor ^ (Ao/a)"^. 

To summarize this analysis: On the one hand we expect 
— ric oc a/Ao, characteristic of the critical behaviour in 
a uniform system. On the other hand — iVc oc (a/Ao)^ 
is dominated by the effect of the jic shift on the density out- 
side the critical region. The latter result was observed in |27|; 
the former cannot be experimentally verified without a direct 
probe of the local density in a 3D cloud. 

By studying the condensed fraction /o at the MF-predicted 
critical point we overcome the problem of the absence of the 
local density probe, and gain more direct insight into the crit- 
ical behaviour in the centre of the trap. Simply put, instead 
of asking how Nc is reduced with respect to by criti- 

cal correlations, we ask how many atoms pile up in the con- 
densate if (at constant T) we increase the total atom num- 
ber to > Nc- Experimentally, the obvious advantage is 
that the condensed and thermal component can be clearly dis- 
tinguished in standard time-of-flight (TOF) expansion, thus 
allowing us to use a "global" measurement technique to ac- 
cess the local behaviour of the gas within the critical re- 
gion. Theoretically, the analogous quantity for a uniform gas, 
no/n (where uq is the condensate density), was first con- 
sidered by Holzmann and Baym OTl . Although the formal 
proof is rather involved, the main scaling result is intuitive, 
no/n (X Auc oc a/Ap |32|. From this result we immediately 
obtain /o oc (a/Ao)^, as illustrated in Fig.[TJb). The presence 
of the harmonic trapping potential still affects the scaling of 
/o with a/Ao, but in this case the results for a harmonic and 
a uniform system are trivially related by the volume of the 
critical region, oc (a/Ao)^. 

To experimentally measure /o we use an optically trapped 
cloud of "^^K atoms in the jT, mi?) = |1, 1) hyperfine state, 
in which the strength of interactions can be tuned via a Fes- 
hbach resonance centred at 402.5 G [331. Our experimental 
system and the procedure for making precise and accurate 
measurements close to the critical point are described in detail 
in ll27ll34ll35 J . Briefly, we prepare partially condensed clouds 
at various values of the scattering length a, and then let the 
number of atoms in the trap gradually decay through inelas- 
tic processes, while finite trap depth and sufficiently high rate 
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of elastic collisions ensure that the sample remains in equilib- 
rium at an approximately constant temperature. For the mea- 
surements presented here, iV w (4 — 5) x 10^, the geomet- 
ric mean of the trapping frequencies in our nearly isotropic 
trap is (I;/27r « 80 Hz, and T « 250 nK, corresponding to 
Ao ~ 10"* aQ, where ao is the Bohr radius. 

To discern condensed fractions as low as ~ 0.1 % we turn 
off the interactions (switch a close to zero) during TOF, thus 
minimizing the condensate expansion ll27ll . We numerically 
calculate using standard MF theory (see also |30|) and 
measure the condensed atom number A^o at the point where 
the total atom number is TV = . To eliminate various 
sources of a-independent systematic errors (including abso- 
lute N and Cj calibration) we perform "reference" measure- 
ments in a weakly interacting gas with a/Ao ~ 0.005 |27|. At 
this reference point the expected value of /o is < 10^^ (see 
below), and we neglect it in our analysis. 

Our experimental results are summarized in Fig. [2] Starting 
at zero for small a (in agreement with MF theory), the con- 
densed fraction /o grows to ^ 1% at a w 350 aq. The use of 
a Feshbach resonance in principle allows us to increase a fur- 
ther, but in the more strongly interacting gases the unfavorable 
ratio of the three-body loss rate to the two-body elastic colli- 
sion rate precludes reliable equilibrium measurements |27|. 

We fit our /□ data with a function {a/X^Y where a; is a free 
parameter The fit yields x — 3.9 ± 0.4, in agreement with the 
predicted x = 4. This confirmation of the expected scaling of 
/o with a/Ao is the first main result of this paper 

We now quantitatively relate our measurements to Monte- 
Carlo (MC) calculations for a uniform gas. Following 1281 we 
first define the reduced chemical potential 



X = 



(4) 



327r3(a/Ao)2fcsT ' 
Next, following 1 1 1 1 we calculate Xq, the value of X in the 



centre of the trap for N = (due to logarithmic correc- 

tions this is a slightly different condition from /i(0) = , 
but this distinction is not experimentally observable). We use 
the experimental value 62 = 42 ± 2, and = 11.7 ± 0.1 
GTllMlto get 



3C(3) 
327r3C(2) 



(62 - h 



= 0.067 ±0.005 



(5) 



For a uniform system the reduced condensate density f{X), 
defined by n^Xl = 167r3(a/Ao)/(X), was tabulated in ||23 
using MC simulations. Invoking LDA, for a harmonically 
trapped gas we get 



iVo 



^/2{'i■KY 



Xo 



f{X)^Xo-XdX . (6) 



Writing {Nq/N^Y^^ = oei^'/Xo) and numerically evaluating 
the integral in Eq. (j6]l, using the results of |28|, we get the 
numerical coefficient auc = 10-4 ± 0.4. 

In Eq. ^ No is calculated at N = A^mf normalised 
to A"^. This expression therefore differs from /o by a factor 
N^'^/N°. This difference does not affect the leading (a/Ao)'* 
term and is relevant only at the {a/XoY level. Nevertheless, 
for a direct quantitative comparison, in Fig. [3] we normalise 
the measured Aq values to A"^, and assume the quartic de- 
pendence on a/Ao. The linear fit to {Nq/N'^Y^^ yields the 
experimental value ctoxp = 10.3 ±0.3, in excellent agreement 
with the Monte-Carlo result. 

For another comparison, it is interesting to convert X^ into 
A'o using the standard Thomas-Fermi (TF) law. This MF law 
is valid well below Tc, where A^o ~ but should not hold 
close to the critical point. For a given Xo, the TF law also 
predicts A'o oc (a/Ao)^. However it corresponds to f{X) ~ 
X and gives q;tf = 8.2 ± 0.4. This result underestimates 
the condensed fraction /o by a factor (q:mc/Q!tf)'' ~ 2.6, 




FIG. 2. (color online) Condensed fraction of an atomic gas induced 
by critical correlations. The condensed fraction /o is measured at the 
point where the total atom number is A = A^''^ > Ac. A fit to the 
data (solid line) with the function /o oc {a/XoY gives an exponent 
X — 3.9 ± 0.4, in agreement with the predicted a; = 4. Vertical 
error bars are statistical and horizontal error bars reflect the 0.1 G 
uncertainty in the position of the Feshbach resonance. The insets 
show representative column density distributions after 19 ms TOF. 



FIG. 3. (color online) Comparison with Monte-Carlo calculations for 
a uniform system. To quantitatively compare our data with the MC 
simulations we plot {Nq/N'^Y^'^ versus a/Ao (see text). A linear fit 
gives a gradient of Ocxp = 10.3 ± 0.3, in excellent agreement with 
the prediction Qmc ~ 10.4 ± 0.4. The error bars are obtained using 
the limiting values from Fig. |2j the points with large error bars do not 
significantly affect the fitted value of Qcxp but are clearly consistent 
with it. 
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and we experimentally exclude it at about 4 sigma level. This 
confirms that near Tc mean-field theory fails on both sides of 
the critical point. 

In conclusion, we have studied the condensed fraction of 
an atomic gas induced by inter-particle correlations at a point 
where no condensate is predicted by mean-field theory. Build- 
ing on the recent observation of correlation effects on the con- 
densation temperature of a trapped gas, this work makes a 
more direct connection with the critical behaviour in a ho- 
mogeneous system. We experimentally confirm the predicted 
scaling /o oc (o/Aq)^, which highlights the conceptual dif- 
ference between the interaction shifts of the critical density 
(characteristic of a uniform system) and the critical atom num- 
ber in a harmonically confined cloud. Moreover, we demon- 
strate excellent quantitative agreement between our experi- 
ments and Monte-Carlo simulations for a homogeneous gas. 
In a more general context, this provides an example of the po- 
tential of ultracold atomic gases for quantitative quantum sim- 
ulation of intricate beyond-mean-field phenomena in uniform 
many-body systems. 

We thank J. Dalibard for comments on the manuscript. This 
work was supported by EPSRC (Grant No. EP/1010580/1). 
R.P.S. acknowledges support from the Newton Trust. 



[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys 80, 885 
(2008). 

[2] T.-L. Ho and Q. Zhou, Nature Physics 6, 131 (2010). 

[3] S. Nascimbene etai. Nature 463, 1057 (2010). 

[4] T. Yefsah etal., arXiv: 1106.0188 (2011). 

[5] T. D. Lee and C. N. Yang, Phys. Rev. 105, 1 1 19 (1957). 

[6] T. D. Lee and C. N. Yang, Phys. Rev. 112, 1419 (1958). 

[7] M. Bijlsma and H. T. C. Stoof, Phys. Rev. A 54, 5085 (1996). 

[8] G. Baym etai, Phys. Rev. Lett. 83, 1703 (1999). 

[9] M. Holzmann and W. Krauth, Phys. Rev. Lett. 83, 2687 (1999). 
[10] J. D. Reppy et ah, Phys. Rev. Lett. 84, 2060 (2000). 
[11] P Arnold and G. Moore, Phys. Rev. Lett. 87, 120401 (2001). 
[12] V. A. Kashurnikov, N. V. Prokof'ev, and B. V. Svistunov, Phys. 
Rev. Lett. 87, 120402 (2001). 



[13] M. Holzmann, G. Baym, J.-P. Blaizot, and F. Laloe, Phys. Rev. 

Lett. 87, 120403 (2001). 
[14] G. Baym etal, Eur. Phys. J. B 24, 107 (2001). 
[15] H. Kleinert, Mod. Phys. Lett. B 17, 1011 (2003). 
[16] J. O. Andersen, Rev. Mod. Phys. 76, 599 (2004). 
[17] M. Holzmann etal, C. R. Physique 5, 21 (2004). 
[18] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Phys. Rev. A 54, 

R4633 (1996). 
[19] J. R. Ensher et ai, Phys. Rev. Lett. 77, 4984 (1996). 
[20] M. Houbiers, H. T. C. Stoof, and E. A. Cornell, Phys. Rev. A 

56, 2041 (1997). 

[21] M. Holzmann, W. Krauth, and M. Naraschewski, Phys. Rev. A 

59, 2956 (1999). 
[22] P Arnold and B. Tomasik, Phys. Rev. A 64, 053609 (2001). 
[23] R Gerbier etai, Phys. Rev. Lett. 92, 030405 (2004). 
[24] M. J. Davis and P B. Blakie, Phys. Rev. Lett. 96, 060404 

(2006). 

[25] O. Zobay, Laser Physics 19, 700 (2009). 

[26] R. Meppelink et al, Phys. Rev. A 81, 053632 (2010). 

[27] R. P Smith, R. L. D. Campbell, N. Tammuz, and Z. Hadzibabic, 
Phys. Rev. Lett. 106, 250403 (2011). 

[28] N. Prokof'ev, O. Ruebenacker, and B. Svistunov, Phys. Rev. A 
69, 053625 (2004). 

[29] Note that B2 is not just a constant but contains logarithmic cor- 
rections in a/Ao 1 1 1 1. We neglect these in our discussion since 
they are not discernible at the current level of experimental pre- 
cision. 

[30] For quantitative comparison of theory and experiment we take 
into account that: (i) 62 includes a small MF contribution, 
^MF _ wj -i^Q 1^ which we obtain by a direct self-consistent 
numerical simulation, and (ii) By fitting the data of |27 | with 
the linear shift fixed to 61 = —3.426 we improve the exper- 
imental estimate for the quadratic shift to 62 = 42 ± 2. These 
coiTections do not affect our general scaling arguments. 

[31] M. Holzmann and G. Baym, Phys. Rev. Lett. 90, 040402 
(2003). 

[32] Strictly speaking this linear scaling holds for any fixed distance 
from the critical point given by — fic) (Ao /a)^ = const L31J . 
By applying it to the MF critical point we again neglect the the- 
oretically expected logarithmic corrections to fi^^ — fic, which 
are so far not experimentally observable. 

[33] M. Zaccanti et ai. Nature Physics 5, 586 (2009). 

[34] R. L. D. Campbell et ai, Phys. Rev. A 82, 063611 (2010). 

[35] N. Tammuz etai, Phys. Rev. Lett. 106, 230401 (2011). 



